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Abstract 

We have studied the weak decay constants and the Isgur-Wise form factor of 
the B and D mesons in the heavy quark limit, by employing a relativistic quark 
confinement model. It is an attempt to improve our previous work within the same 
line of thinking, but by incorporating a couple of novel aspects. First, the infrared 
behavior of the heavy quark is considered by modifying its conventional propagator 
in terms of a single parameter v. Second, the mass difference of the heavy meson 
and heavy quark: E = ran — Mq has been included. It is found that the weak 
decay constants depend strongly on the mass difference E with a relatively mild v 
dependence. As for the Isgur-Wise function it is controlled more sensitively by the 
infrared parameter v, leading to its suppression at maximum meson recoil. 
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1 Introduction 



The theoretical treatment of heavy-quark hadrons: those containing a heavy (b or c) 
quark (Mq >> Aqcd), is considerably simplified due to a new heavy quark (spin-flavor) 
symmetry Q. Under this symmetry a heavy-quark meson may be pictured as a system 
of a cloud of light degrees of freedom (light quarks and gluons) surrounding the almost 
static color source (a heavy quark) whose four-momentum being quite close to its on- 
mass-shell value. In practice, it allows for expressing all the B — > D(D*)lu decay form 
factors in terms of a single universal form factor: the Isgur-Wise function (IWF), which 
is determined primarily by the dynamics of the light degrees of freedom. 

The systematical study of this symmetry is provided by the heavy-quark effective 
theory (HQET) developed by Georgi et al. 0. In HQET the heavy degrees of freedom are 
integrated out, while their effects as virtual particles are taken into account by introducing 
non-renormalizable higher- dimension operators proportional to the inverse powers of the 
heavy quark mass ||. However, IWF cannot be directly obtained within this scheme 
since it is essentially determined by the physics in the non-perturbative domain. Thus, 
to work this quantity out various quark models M , Q , M , [0] , M , [0,0], QCD sum rules 
I2| , , |TJJ , |T5j , and lattice calculations |I7| have been employed so far. The present 



work is devoted to the extension of our previous calculation of IWF within a relativistic 
confined quark model [0. To constrain the model parameters we will also calculate the 
weak decay constants of the B and D mesons. In order to set a basic stage for our study, 
we shall first briefly review some of the model results to date in the following. 

In the nonreltivistic confined quark model j|] IWF is expressed in terms of the overlap 
integral of the initial and final meson wave functions of the Gaussian form which are 
determined by fit to meson spectra (note that heavy quark mass is assumed to be infinity 
here). With w = v ■ v' where v(v') is the four- velocity of the initial (final) heavy-quark 
meson, it reads 

£H = exp{-pV-l)}- (1) 

The slope parameter obtained from the meson spectroscopy is p 2 = — £'(1) = 0.33. An- 
ticipating that the result would be insufficient for large recoil, the authors of have 
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introduced a universal reduction factor (from a fit to the pion electromagnetic form fac- 
tor) which has resulted in an increase in the slope parameter: p 2 « 0.64. We note here 
a recent claim 0] that a due consideration of the relativistic recoil correction and the 
Wigner rotation of the spin of the light quark should increase the slope parameter to 
greater than unity. 

To improve certain aspects of nonrelativistic approachs a simple relativistic-oscillator 
model with the light-front variables was adopted to find IWF ||, 

£H = ^-exp{-(2p 2 - 1)^}. (2) 
w + 1 w + 1 

Here, the slope parameter is related to the transverse momentum distribution of the light 
quark inside the heavy meson. By comparing with the data its value has been extracted 
as p ~ 1.14 ±0.23. 

We next discuss other quark model approaches with or without due consideration to 
confinement. First, a simple model based on quark loop graphs with an ad hoc form factor 
at each heavy-meson heavy-quark light-antiquark vertex was proposed in Ref. ||. 
It was found that the physical quantities (including IWF) are quite sensitive to the mass 
difference; E = itlh — Mq, of the heavy meson and heavy quark when it exceeds m q : the 
light-quark (constituent) mass of the model. Then the mass function (self energy), decay 
form factor (IWF), etc. of the heavy meson made up of a heavy-quark light-antiquark 
loop develop an imaginary part which becomes of the same order of magnitude as the 
real part. This is of course due to the absence of the quark confinement in the model. 
Eventually the imaginary part has been discarded completely by arguing that its neglect 
would mimic the quark confinement. To date, the justification of such a procedure is not 
known. 

A modified MIT bag model was applied to the calculation of IWF for the ground state 
to ground state semi-leptonic decays involving the b — > c transition |L0|1 . The evaluation 
of the overlap integral yielded a result, which is claimed to be well approximated by the 
following form |13 
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The heavy-quark meson wave function was obtained by solving the bound state Bethe- 
Salpeter equation with a modified one-gluon-exchange ladder approximation JO]. The 
decay constant fs and IWF were then calculated from the wave function. The slope 
parameter at zero recoil was found to be rather large: p 2 = 1.8 — 2.0, whereas the weak 
decay constant turned out very small, ~ 50 MeV. Presumably, the one-gluon-exchange 
may not be very realistic for generation of the heavy meson bound state (this interaction 
alone does not confine the quarks). 

Our previous approach for IWF was based upon the quark confinement model 
(QCM) which incorporate the confinement of light quarks by devising a quark propagator 
that has no singularities, thus forbidding the production of a free quark [^5j . This model 
then allows one to perform covariant calculations of Feynman diagrams with dressed light 
quark propagators free of simple pole. The calculation has given, in the heavy quark 
limit: Mq = mjj — > oo, the IWF to decrease slower than that obtained by many other 



approaches except the bounds []I9| to be discussed below, 

1 



l + R 



2 1 , 

R + — == h\[w + Vw 2 - 1] 



(4) 



1 + W y/w 2 - 1 

with the parameter fixed by the light quark confinement function (see the next section) to 
be R = 0.65. The corresponding slope parameter was extracted to be equal to p 2 = 0.43, 
quite small as compared with the valued obtained from a fit to the data by the relativistic 
oscillater model as discussed above. It should be mentioned that the IWF was obtained 
also with finite heavy quark masses, and the result turned out to be very close to the one 
in the heavy quark limit. 

It should be useful at this point to refer to a rather peculiar model [[J which, in part, 
has motivated our present investigation. It is based upon the observation of the possible 
infrared behavior of the heavy quark propagator. The physics here is not quite dynamical 
but mostly geometrical, i.e. in terms infraparticle (heavy quark) propagation only without 
taking into account the light degrees of freedom. The obtained IWF has the form 

£0) = l-exp {-—== , (5) 



Vw 2 - 1 . 

with the slope at zero recoil equal to zero. The model does not appear particularly 
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realistic, but its consideration of the infrared propagator is quite novel which we have 
adopted in our present work. 

Different from various quark model approaches discused so far, the QCD (spectral) sum 
rule [0 appears to be the only alternative so far in effectively dealing the long distance 
physics of the heavy mesons, aside from the direct numerical simulation by lattice QCD. 
This approach is based upon the hypothesis of quark-hadron duality in the calculation of 
physical correlators and is combined with HQET to simplify the physics. After several 
simplifying assumptions the first sum rule calculation ]12[ found 



£0) « exp{-0.37Vw 2 - 1} (6) 

with the infinite slope parameter. 

Later, this approach was improved |TJ| ||15| , by including the next to leading 



order renormalization group improvement, finite heavy-quark mass effects, etc. It was 
found that these effects are quite large for the weak decay constant and makes the D- 
meson as much as about 40%, |13| causing a large deviation from the asymptotic scaling 
law foy/mD = fBy/mB- Also a rather strong dependence of this quantity on the mass 
difference E = itlh — Mq was observed . As for IWF the E dependence was fairly weak or 
cancel out completely. The main uncertainty brought in the calculation of this quantity 
comes from the QCD perturbative estimate of the higher-resonance contributions to the 
spectral function in the dispersion integral representation of the hadronic correlator. It 
turned out that the different choice of the spectral function and of its integration domain 
produces the values of the IWF form factor at maximum recoil varing in the region 
0.40 < £(w max ) < 0.75. It was also shown that the infinite slope parameter of [JT2J] has 



resulted from a naive choice of the integration domain of the spectral function. There 
appears, however, no single best choice for this domain. The "best" result compared with 
the experimental data JJ7) may be parametrized as in eq.(3), giving p 2 = 1.3. 

The lattice QCD has been exploited by two groups |L6[ and flTfl in the study of 



heavy meson physics. The IWF was obtained by calculating the elastic D amplitude (or 
form factor) < D'\c r y^c\D >, in the quenched approximation. Fitting the result to the 
relativistic-oscillator parametrization Eq.(2) gives 
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1.41 ±0.19 ±0.41 @ 



P 



1.211 
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Lastly, we should mention that the bounds on the IWF were studied without relying on 
particular models, but only with the consideration of analyticity and dispersion relation 
for the heavy meson elastic form factor assuming the heavy quark symmetry and the 
hardon-quark duality It was claimed that the bounds is very strongly constrained. 
These bounds prohibit the IWF from falling by less than 12% and more than 21% from 
zero to maximum recoil. The lower bound conflicts with many model predictions. It has 



been argued by |£0| , |2l] , |22| , and |23) that the careful inclusion of heavy quark- ant iquark 
bound states, lying below the heavy meson pair production threshold, may gives no such 
small lower bound. The authors of |0| recently included the effect of the T states below 
the BB threshold claiming that the rigorous upper bound for the IWF slope at zero 
recoil is now p 2 ~ 6.0, although, on phenomenological grounds, it may be diminished to 
P 2 < 1.7. 

The main goal of the present paper is to improve our previous work [[7| as we believe 
that our QCM approach possesses a sound dynamical background regarding its modelling 
of confined light quarks, successful in many static and non-static properties of the light 



hadrons |25|. Thus the present investigation will focus on how one should treat the 
behavior of the heavy quark within the context of this QCM. 



2 Quark propagators 

The quark confinement model (QCM) |25[ is specified by the interaction Lagrangian 
inferred from the partial hadronization of the DCQ action. Then the transition between 
heavy meson H and heavy Q and light q quarks is described by 

£h(x) = g H H{x)Q{x)T H q(x) (7) 

with the coupling constants gu defined by what is usually called the compositeness condition 
which means that the renormalization constant of the meson field is equal to zero: 
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Z H = l-g 2 H Tl' H (m%) = 0. (8) 

Here Tl' H is the derivative of the meson mass operator. The compositeness condition also 
provides the right normalization of the charge form factor Fjj(0) = 1. This could be 
readily seen from the Ward identity 

AW) = ti^^f 1 = 9 2 H^TUP,P) = F H (0) = 1. (9) 
with the three-point function T M (p,p') defined as 

T£ H (P,P') = 4^/ ^{ VS Q^- t)-fS<iW- ^S q m). (10) 
The matrix elements describing the leptonic decays H(H*) — > ev (H = B,D), and 
semileptonic ones B — > D(D*)eu are written down as 



where 



and 



where 



M"{H(H*)^ev} = g H T£(p), (11) 



TM = h J Bi • tr {°" w ~ (12) 



M B^D(D*)eu = 9Bg D (D*)T£ H ,(p,p'), (13) 



t£ h ,(p,p') = ^- 2 J ^r{ 7 5 ^(^- f)<rs Q {ft- firsM}. (u) 



Now we discuss the forms of quark propagators. 
(1) Confined light-quark propagator. 
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Here we briefly explain the salient feature of QCM concerning the light-quark propaga- 



tor. The principal assumption |25[] is that the propagator is an entire analytical function 
in the complex-" p 1 " plane with the constituent mass m q = zA being smeared by complex 
measure da, such that 



da, 



Az- 



1 

A 



pf 

A2 



A 



■pf 
A2 



(15) 



with the confinement functions a(u) and b{u) defined by 



a(u) 



da, 



z 2 + u 2 



b(u) 



da, 



z 2 + u 2 



(16) 



The scale parameter A characterizes the size of the confinement region. It has turned 



out |25| that the low-energy physics depends only on those quantities which involve the 
integrals of a(u) and b(u) together with u a (a > 0) but not on the detailed shape of 
these functions. We have devized a simple choice of the confinement functions ||2|| : 



a(u) = a,Q exp( — — a\u) b{u) = 6oexp(— u 2 + b\u). 



(17) 



The parameters a^, bi, and A have been determined from the best model description of 
hadronic properties at low energies and the following values were found: 



a = b = 2, ai = 1, b x = 0.4, 

(18) 

and A = 460MeV, 
which describe various physical observables quite well [^J. 

(2) Infrared heavy-quark propagator. 

Since a heavy quark (with mass Mq) in a heavy meson is under the influence of soft 
gluons (which sets the scale Aqcd « Mq), it may be regarded as nearly on its mass-shell 
where the infrared regime should take place for its propagation. The infrared behavior for 
one-fermion Green's function (propagator) has been investigated in various papers (see, 
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for instance, 0, |26|],and the references therein). The result is well-known only for abelian 
gauge theories: 



S(p) ~ (m- i>-ieY x ~\ (19) 

where v = (as/47r)(3 — A) with A being the gauge parameter. It seems quite resonable to 
see how sensitive the physical observables may be to this modification even though the 
integrals in Eqs.(lO), (12), (14) have no infrared divergences. To make our present study 
simple, we choose the heavy quark propagator as in Eq.(19), 



S Q (k + p) = (M Q - ft- fr 1 -" (20) 

assuming that v is a free parameter. Using the conventional notation for the heavy meson 
four-momentum p = (Mq + E)v, with v being the four-velocity of the heavy meson (or 
of the heavy quark since there is a velocity selection rule: 0, one readily finds that this 
propagator satisfies the form necessary for the correct heavy quark limit: 



S Q (k + p) 



1 



Mn- U- 



1 1+ " + (T L). (2!) 



" 2(-E-kv) 1+u ' ~ y M ( _ 
The physical observables are defined by two parameters: the heavy meson heavy quark 
mass difference E = mn ~ Mq which typifies the binding of the light and heavy quarks, 
and v characterizing the infrared behavior of the heavy quark as disscussed above. It 
should be mentioned here that in our previous work j7| we have not considered the effect 
of the non- vanishing E. In the heavy quark limit both the meson mass mn and heavy 
quark mass Mq are sent to infinity, so it does not appear to make sense to keep this value 
finite. However, it characterizes the binding of the light and heavy quarks, thus an very 
important parameter particularly in constraining the values of the weak decay constant, 
as we shall see later. We also note that this quantity plays a significant role in other 
approaches. 



3 Weak decay constants and IW form factor 

The calculation of the two- and three-point functions (see Eqs.(10),(12), and (14)) is 
considerably simplified by taking only the leading order in the 1/Mq expansion (the 
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heavy quark limit) since the calculations of the traces and the integrations over internal 
momentum k factorize in this limit due to the following identities: 



tr [0 M (1+ flr(z+ #)] = tr [0"(1+ f)F(z+ ^kv)] = tr [0"(1+ f)T] (z - kv) (22) 



tr [T(l+ ^)7 M (1+ tfp(z+ ft)] = 2tAr [rr] (z - kv) 



(23) 



tr [r(l+ /)0"(1+ ^)7 5 (^+ 

tr r(l+ ^)0"(1+ ^) 7 5 ^ + ^ + ^ ( ^ + ^ _ ^zM_^_ ^ 



(24) 



tr [r(l+ /)0"(1+ ^) 7 J 



2(w + l) 



2(w- 1) 



(iu + i) ; 

where w — v • v'. Here we have taken into account that {r, ^} = for both T = 7 5 and 
7", which follows from the transversality of the vector field. 
Using these identities in the Eqs.(lO), (12) and (14) gives 



T HH'(P,P') 



2p 



= A 



Here we use notation: 



\m H J \2A U J 47r 2 



1 ^ 3 Jt\j,l)-L , 



A 



' 3 T^( E \ T» 

; J 2 Ix^'^i' 



m H J 2AMtt 2 z V A 



A 



J. 



(") 



,E 



(2A^) 2 4vr 2 d V A 



,w) ■ L 



2 • 



L = < 



/ (r = 7 5 ) 
I s a/3 (r = r) 



(25) 
(26) 
(27) 



jf (r = 7 5 ) 

^ Me' 1 (r = 7 a ) 
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L 2 



(v + v'Y (r = 7 5 ) 



where e is the polarization vector. The structural integrals and J% are defined as 



# (EM = /^/^_^Z^^^ (29) 

and they are calculated in the Appendix. 

The i^Qg coupling constant is determined from Eq.(8) and may be written down as 



9 (m H ) = 2A^^ 



(30) 



Finally, the matrix elements of the leptonic and semileptonic decays are written as 

M B^D{D*)ev = y/m B m D £{w) ■ L% (32) 
with the weak decay constant /# and the IWF given by 



^ - rjp(E/A, i) 



(33) 



= (34) 



jf } (£/A, W ) 

Now we discuss the dependence of physical observables on parameters E and v. Note 
that hereafter E is given in units of A (= 460 MeV), At this point, it may be useful to 
stress that unlike in Ref. |§ there is no imaginary part in our matrix elements for any 
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value of E since there is no Qq production threshold due to confinement. The value of 
E is not arbitrary, however. The first restriction comes from the requirement that the 
integral J^(E,1) in Eq.(30) should be a positive definite quantity. Since this integral 
for v > is defined in terms of the derivatives of the confinement functions, it develops 
zeros above certain values of E. The behavior of J%(E,1) for v = 0,0.5,1 is shown in 
Fig.l. Another restriction on E comes from the value of the heavy meson weak decay 
constants. From eq.(33) it is clear that our result in the heavy quark limit obeys the 
scaling relation fs\ /m B = fD^/ m D- The dependence of the weak decay constant values 
on E is shown in Fig. 2, and it is clear that this quantity is very sensitive to the change in 
E: strongly increasing function of E. This is in line with what was observed in QCD sum 
rule approaches, see for example, |13[ . In the present result the E=0 case is unrealistically 
small, suggesting that even in the heavy quark limit E must be kept finite as it should be. 

The values of the weak decay constants obtained in different models are in variance 
but eventually f D never exceeds ~ 250 MeV @, J27|, |2S|, @, H when the pion decay 
constant is normalized to f n =132 MeV. We will use the above observations as a guide to 
constrain E. The numerical results for the upper bound obtained in this way for various 
choice of the infrared parameter v are shown in Fig. 3. It is interesting to note the lower 
bound for E obtained in Ref. E> 0.5, which might eventually be implemented to get 
the restriction on the infrared parameter v which turns out as v < 0.75. Note, however, 



that there is an antithesis to the way this lower bound has been derived fl32 |. 

The IWF is plotted in Fig. 4 for various values of E and v. One can see that the E- 
dependence of this function is much weaker than that seen in the weak decay constants. 
This agrees qualitatively with the result obtained by using QCD sum rules fll3f , fl4 |, f!5 |. 



Note, in particular, that in [13[ the IWF is completely independent of E, which is rather 
curiuous intuitively. In the present result the E-dependence decreases for increasing v: 
for v = and at maximum recoil the difference between the values of the IWF for E=0 
and E=1.0 is no greater than 12%. The general tendency is that the IWF is suppressed 
for larger values of E. 

On the other hand, the dependence on v is much stronger as seen in Fig. 5. For 
examlpe, at maximum recoil the IWF is suppressed by 13% for v = 0.5, and by 23% for 
v — 1 in comparison with v = when E=0. A similar tendency is observed for other 
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values of E. 

The E- dependence of the slope parameter at zero recoil: p 2 = — £'(1), which is an 
important quantity characterizing the IWF, is shown in Fig. 6. One can see that it is 
contained within the following region: 



0.42 < p 2 < 0.64 (u = 0), (35) 

0.64 < p 2 < 0.82 [y = 0.5), (36) 

0.62 < p 2 < 0.64 (u = 1). (37) 

The result is consistent with the generous Bjorken lower bound: p 2 > 0.25 [33]: but 



somewhat less than the values obtained in other approaches (see, Table 2). It should be 
noted, however, that a recent fit to the CLEO B — > D*lu data with various analytic forms 
adopted for the IWF give p 2 ~ 1.1 ±0.5 ±0.3, consistent with our present result |53j] with 
v ph 0.5. Also note that this quantity does not necessarily determine the global behavior 
of the IWF, as seen in Fig. 8. Also, its accurate determination from experimental data is 
not trivial, in good part, due to the small count near zero recoil. The effort is being made 
from the CLEO data for the B — > D*W and we await a quick release of the analysis. 

Plot of the IWF for the set of parameters {is, E} providing the reasonable values for 
the weak decay constants (see, Table 1.) is presented in Fig. 7. One can see that, for 
example, the curve with the parameters v = 0.75, E=0.5 (or E=230 MeV) is consistent 
with the experimental result within 19% off the best fit taken from Ref. []38|| . In Fig. 8 
the IWF taken from the other approaches are shown. 

4 The effects of finite masses in the weak decays 

To study the correction to the heavy quark limit, we will take into account the effect 
of finite heavy quark masses in the weak decays of B and D mesons. This is of special 



interest since the QCD sum rule fl3]j, [E^j] and the lattice calculations [E3|, [E9| give an 



indication of a strong breakdown of the scaling law: 
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f D = 170 ± 30 MeV f B = 190 ± 50 MeV 







f D = 150 -> 195 MeV / B = 150 — > 185 MeV 



/ D = 200 ± 25 MeV / fl = 230 ± 10 MeV 



28 



/ D = 185ifjf MeV /, 



160±f±f| MeV 



/d = 208(9) ± 35 ± 12 MeV / B = 187(10) ± 34 ± 15 MeV 



ISO 



In [27]] the spread in the value of the decay constants originates from varying the values 
of the heavy quark masses. Also the results are found without adopting any HQET, and 
the sum rule result with the method of Hilbert moments not listed here has given somewaht 



different values. Note that in |13| and |28[ the complete break down of the scaling law 
even brings the value of fo lower than that of 

To simplified the calculations, we consider only the case with v = in view of its very 
mild influence on the decay constants discussed earlier. Then we have 



^ ^4 0) (M Q /A,E/A) 



(38) 



with the structural integrals defined as 



4°\M,E) 

4°\m,e) 

where 



dtt 



1 + t/Mf 



t(l + J rT )-E{l+ h 



2M' 



2M' 



, . 1 * 



2 1 + t/M 



b(z) 



(39) 



dtt r . . t 

; — —< a(z) H — 

;i + t/M) 3 l K J 2 



3 + — 



t (1 + E/Mf 



M l + t/M 



b(z) 



l + t/M 



t-2E(l + — )}. 
v 2M>) 
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It is easy to show that eq.(33), whose structural integrals being defined in the Appendix, 
is reproduced in the heavy quark limit Mq — > oo. 

We have adopted the values of the heavy mesons: m D = 1.87 GeV, m B = 5.27 GeV, 
and varied E = ran — Mq which, in turn, determines the finite heavy quark masses. Just 
like in the heavy quark limit as discussed earlier, /# has turned out quite small for small 
E, and the effect of the finite Mq is found rather insignificant. However, for E > 0.5 f B 
increases significantly. In order to compared with the result in the heavy quark limit as 
listed in Table 1, we we present here the finite quark mass result for E=l: 



f D = 175 MeV f B = 125 MeV, (40) 

which is quite reasonable as compared with the lattice and QCD sum rule results quated 
above except that the B-meson decay constant is somewhat on the smaller side. When 
compared with our value in the heavy quark limit, the reduction is about 10% for the B 
meson and as large as 25% for the D meson. So obviously this effect is important and 
breaks the scaling law. The percentage reduction here is in agreement with what was 
found in ref. [13] within the QCD sum rule when the extra reduction factor (Mq/itlh) 1 ' 5 
in this work is left out. We should note that it is this extra factor in this reference 
that eventually reverses the magnitudes of /o and f B (note that there the adopted value 
of E is somewhat smaller for D), thus consistent with the lattice result in |[28|| . In our 
present approach this reduction factor is apparently absent: we wonder if it is a model 
independent factor we have missed, or it is due to the specific way the parameters were 
defined in performing the Laplace sum rule in the heavy quark limit in |TB[. In any event, 
we should stress that the present section is not intended as predicting the values of the 
heavy meson decay constants: rather it is to demonstrate the importance of the finite 
heavy quark mass effect to correct the heavy quark limit as long as this specific physical 
quantity is concerned, particularly for the D meson. 

As for the finite quark mass effect on the IWF, we already studied it in our previous 
work @ where E = 0. The effect was found to be rather small and we expect a similar 
conclusion to hold here, thus has not been investigated. 
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5 Conclusion 



As an extension of our previous work |7| , we have presented a quark model approach to the 
decay of heavy quark mesons where a confined light quark and the near on-shell infrared 
behavior of the heavy quarkpropagation are taken into account. The infrared behavior of 
the heavy quark is modelled in terms of a single parameter v which modifies the simple 
free Feynman propagator. The study has been carried out in the heavy quark limit: 
Mq — > oo. However, we have retained the finite value for E = ran — Mq: the difference 
between the heavy meson and heavy quark masses. We have found that the weak decay 
constant of the heavy quark mesons depends strongly on E while its //-dependence is 
rather mild. This fact has then been used to constrain the value of E by imposing a 
conservative upper bound for fn to be ~ 250 MeV for the D meson decay. In this way 
it is found that for a reasonable value E ~ 460 MeV with v =0 ( corresponding to a 
free heavy quark propagation) the value fo is found to be consistent with approaches like 
QCD sum rule @, j27| and a lattice simulation j29|, g. 



In the same context we have studied the universal semi-leptonic decay form factor: the 
Isgur-Wise function. Unlike in the case of the decay constants, it is found to be controlled 
more dominantly by the infrared parameter v which gives a desirable suppression of this 
function for v > 0. As an example, for a set of parameters: E=0 and v = 0.75 giving 
reasonable values for the weak decay constants, the I WF turns out to be consistent with 
the best empirical fit to the ARGUS data |37| nearly within the experimental error bars: 
we may need a little extra suppresion. 

Our treatment of the infrared behavior of the heavy quark propagator in terms of a 
single parameter may seem too simple. It is motivated by studies in the abelian gauge 
field, thus is not rooted from real QCD. However, we believe that it certainly points to 
the direction of further investigation, as it appears reasonable enough to believe that a 
free heavy quark propagator (with a simple pole at the quark mass) may not be realistic 
but needs to be modified in such a way to respect the heavy quark symmetry. 

We have also studied the effect of the finite heavy quark mass, and found that it gives 
a significant deviation from the scaling law for the weak decay constants obtained in the 
heavy quark limit (and also in the non-relativistic quark models). We thus confirmed the 
claims from certain lattice and QCD sum rule approaches. 
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To conclude, we want to reiterate the following: while the importance of retaining the 
finite value for E appears to be established in other approaches, the possible significance 
of the infrared behavior of the heavy quark propagator has not been paid much attention 
except in a simple model II. While we ourselves are in persuit of a more realistic form 
of the heavy quark infrared propagator, we will certainly feel happy if the present work 
may invoke certain interest in this direction. 
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APPENDIX 



Two-point function. 
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Three-point function. 
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with W given by 
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(l)i/ = 



W = l + 2r(l -r)(w- 1). 



1 CO 

dr 
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o o 
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v = 0.00 
E=460 MeV 


v = 0.50 
E=276 MeV 


v = 0.75 
E=230 MeV 


f D (MeV) 
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216 


226 


fa (MeV) 


137 


140 


135 


P 2 


0.61 


0.83 


0.92 



Table 2. 



Bjorken [|33j 


>\ 


Isgur § 


0.63(0.33) 


Rosner 


35] 


1.44 ±0.41 


Mannel 


36 


1.77 ±0.74 


Neubert 
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1.28 ±0.25 


Bernard [|16| 


1.41 ±0.19 ±0.41 


UKQCD Collab. [0] 


1.213 
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Sadzikowski 10] 


1.24 


Kugo 11] 
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